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Abstract 

A basic problem in linear particle optics is to find a symplectic transformation 
that brings a symmetric matrix S (the beam or bunch matrix) to a special 
diagonal form, called normal form. The conventional way to do this involves an 
eigenvalue-decomposition of a matrix related to U, and may be applied to the 
case of 1, 2 or 3 particle degrees of freedom. For 2 degrees of freedom, a different 
normalization method involving ’’real Dirac matrices” has recently been proposed 
m, m- In the present article, the mathematics of real Dirac matrices is presented 
differently. Another normalization recipe is given, and more general decoupling 
problems are solved. A 3D visual representation of S is provided. The 
corresponding normalization method for 1 degree of freedom involving ’’real Pauli 
matrices” is also given. 
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1 Introduction 


In accelerator physics, particle beams are usually chopped into bunches of particles. 

The single particle has 3 degrees of freedom (DFs), and occupies a point in 
6 -dimensional phase space. The distribution of bunch particles in phase space can be 
roughly described by a symmetric 6 x 6 matrix S of 2 nd moments, called bunch matrix. 
The evolution of the bunch matrix along a beamline can be described with a symplectic 
6 x 6 matrix R, called transport matrix. This article concerns the basic problem of 
hnding a transport matrix N that will bring the bunch matrix to a special diagonal 
form, called normal form. The conventional way to do this involves an 
eigenvalue-decomposition of the bunch matrix times another 6 x 6 symmetric matrix 7 , 
called the symplectic form [3]. This normalization method is described in section 2. 

Often, only the two transversal particle DFs are of interest. The bunch matrix then 
reduces to order 4, and is called a beam matrix. The normalization problem becomes 
accordingly simpler. For this case, Ch. Baumgarten has recently proposed a different 
normalization method involving ’’real Dirac matrices” n. 0 - In section 5, 
Baumgarten’s method is presented in a mathematically different way: a 2-index 
numbering of the real Dirac matrices is used; more general elementary transformations 
are dehned; the beam matrix instead of a force matrix is transformed; more general 
decoupling problems are solved; and an alternative normalization recipe is proposed. 
Furthermore, general determinant and inversion formulas for square matrices of order 4 
are given, and a visual 3D representation of the beam matrix is provided. 

If only one particle DF is considered, the beam matrix further reduces to order 2. The 
normalization method corresponding to Baumgarten’s method then involves ’’real Pauli 
matrices”. This method is hardly ever used, since the normalizing transformation N 
can be directly expressed in terms of the components of the 2 x 2 beam matrix S [ 1 ], 
[5]. However, getting familiar with it is a good preparation for the more difficult case of 
2 DFs. Normalizing a beam matrix for 1 DF with real Pauli matrices is therefore 
presented in section 4- 

Both the real Pauli and Dirac matrices are real representations of Clifford-Algebras [7j, 
0.0- The details are explained in section 3. 

Introducing real Pauli or Dirac matrices into particle optics means, expressing the 
transformation and beam matrices in a non-standard basis. Since matrix elements in the 
standard basis come with dehnite units (usually mm“mrad^, with a, 6 = ... — 1 , 0 , 1 ...), 
the matrix components in the new basis have mixed units, and have no obvious physical 
signihcance. It’s a case of adding oranges and apples. This means that the use of real 
Pauli or Dirac matrices allows mathematical tricks, but doesn’t add to physical insight. 
Another short-coming of real Pauli or Dirac matrices is, that they cannot be used to 
describe phase vectors, or more generally (phase) tensors of degree other than 2 . 

Finally, there is no matrix algebra of order 6 corresponding to the real Pauli or Dirac 
matrices, since real irreducible representations of Clifford algebras have order 2,4, 8 .... 
In other words, there is no obvious non-standard matrix basis for 3-DF-particle optics. 

The hgures in this article were created with the software package GNU Octave (version 
3.6.4 for the operating system Windows) [B]. I thank the Octave development 
community for making this software freely available, and for their excellent work. 
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2 Linear optics for particles with three degrees of freedom 

2.1 The particle bunch 


A particle bunch can be described by the distribution of its particles in phase space. A 
single particle has 3 degrees of freedom, and the particle phase space is 6-dimensional. 
The phase space coordinate system is chosen co-moving with the bunch, and the 
coordinates and momenta are denoted by Xt = {x,x',y,y',l,6) [ 1 ], [ 5 ]. 

The main information about the bunch is contained in the lowest moments of the 
particle distribution in phase space: the vector X of first moments, describing the 
bunch position, and the matrix S = (X — X)(Xt — Xt) of second moments, describing 
the bunch shape. 


2.2 Symplectic transport 

In this section the following extended versions of the three real Pauli matrices ai 
and ia2 are used: 
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( 2 . 2 ; 


a and (3 are bireal units, 7 is a complex unit, see section 3 . 7 is the symplectic form on 
the phase space. Now let us return to linear particle optics. 

As the bunch runs down the beamline, its particles describe trajectories X in phase 
space, and the bunch mean position X and 2 nd moments S change. The evolution of 
all these quantities is described with a 6 x 6 transport matrix R: 


X = RX( 0 ) 


X = RX( 0 ) 


S = RS( 0 )RT 


( 2 . 3 ) 


The particle motion is such that the phase space volume enclosed by many 
neighbouring particles is preserved, and the transport matrix R is symplectic: 


R7R^ = 7 


or 


R^7R = 7 


or 


R ^ = — 7 R ^7 


( 2 . 4 ) 


Now let r be the logarithm of the transformation R: 

r = In R R = exp r 


(2.6) 


It can be shown that the above condition on R translates to 


r = 7r^7 or 7r = —r^7 or r7 = —7r^ 


( 2 . 6 ) 
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I will call matrices with this kind of symmetry --y-symmetric. The second and third eqn. 
fl2.6p say that a 7-symmetric matrix left- or right-multiplied with 7 becomes symmetric. 
Conversely, a symmetric matrix left- or right-multiplied with 7 becomes 7-symmetric. 
Therefore, to construct a symplectic transformation R, take any symmetric 6x6 
matrix, left- or right-multiply it with 7, and exponentiate. 


2.3 Invariants of the bunch matrix 

Given a 6 x 6-matrix of 2 nd bunch moments S, three numbers can be calculated that 
are invariant under symplectic transport. These are the bunch emittances 61,62,63. 

They are a measure of the bunch size in phase space. 

To hnd the emittances, note that S7 transforms according to 

S7 = RS( 0 ) 7 R-^ ( 2 . 7 ) 

This similarity transformation leaves the characteristic polynomial of 117 invariant. The 
latter is an even function of its argument, as the following calculation shows: 

det (S7 — Al) det (S -|- A7) = det (S — A7) = det (S7 -|- Al) 

The squared emittances now appear as coefficients in the characteristic polynomial: 

det (S7 - Al) = (A^ + e]) (A^ + e]^) (A^ + e]jj) ( 2 . 8 ) 


The eigenvalues of S7 are therefore A = ±*6/, ±ieni- 


2.4 Normalizing the bunch matrix 


[ 3 ] We are looking for a symplectic transformation N that puts the bunch matrix S into 
normal form: 


S = NSNT 
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To this end, right-multiply eqn. fl 2 . 9 p with 7 and use eqn. 02 . 41 ) . The result is 
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This looks similar to the eigenvalue-decomposition 


ie/ 0 0 0 0 ^ \ 

0 ze/ 0 0 0 0 

0 0 -ien 0 0 0 . . 

0 0 0 ien 0 0 ^ ^ 

0 0 0 0 — i^iii 0 

0 0 0 0 0 i^iii j 

where we can normalize the complex matrix E to symplectic. 

The RHS’s of eqn.s fl2.10p and fl2.1ip can be made to agree by means of an intermediate 
symplectic transformation Q: 

N = EQ S7 = Q-^AQ ( 2 . 12 ) 


/ 


E7 = EAE^ 


A = 


A calculation shows that the following Q does the trick: 


/ 1 z 0 0 0 0 \ 

z 1 0 0 0 0 

l + il 3 _ I 001 z 00 
^ “ V2 0 0 z 1 0 0 

0 0 0 0 1 z 

\ 0 0 0 0 z 1 / 


( 2 . 13 ) 


For brevity of notation, let us dehne the matrices with r = I, II or 

III. They are the restrictions of these 6x6 matrices to the 2 x 2 block, e. g. 
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etc. 


To obtain the normal decomposition S = NSNt of a given bunch matrix S, proceed 
as follows: 

1 . Calculate the eigenvalue-decomposition S7 = EAE^^. 

2 . Order the eigenvectors in E such that A = —ia^ej — ia^^en — ioL^^^em. 

3 . Normalize the eigenvectors in E to symplectic: ET7E = 7. 

4 . Set N = E (1 + z/ 3 )/\/ 2 , and S = IR7 + + l^^^eni- 

In step 3 , the real transpose, not the conjugate transpose, must be used. 

In step 4 , N should come out real. 

In appendix A I give a short Octave [6] test program (it should run on Matlab, too) 
that will apply these steps to a random beam matrix. 
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2.5 Invariance group of the bunch matrix 

Which transformations I leave the bnnch matrix S invariant; ISIt = S ? 

If the bnnch matrix is normalized, the invariance transformations are phase rotations 
in the 3 snbspaces {x,x'), {y,y'), 

I = exp 

with ip\ip ^^arbitrary 

If the bnnch matrix is not normalized (which is the general case), let N be a normalizing 
transformation. The gronp of invariance transformations is then given by 

I ('^^, ip^\ = N exp (7^'?/’^ + 
with ip\ip ^^arbitrary 
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3 Matrices and Clifford algebras 

3.1 Definitions 


The real d x d matrices form a non-commutative ring: you can add and multiply the 
matrices. Nowadays mathematicians prefer to talk of an associative algebra: the 
additional property is scalar multiplication, which is the same as multiplication with a 
multiple of the unit matrix. 

A d X d matrix Z = {Z^n) can be written in vector notation as 

Z = + ... + Zd,d-i + ( 3 . 1 ) 

Here, the basis of the d x d matrix algebra is Em,n = dk^mdi,n- 

According to the problem at hand, there may be a better choice of the basis, one that 
simplihes calculation. Particularly attractive are basis matrices that square to ± 1 . The 
general matrix then has the form 

Z = IZq + f3iZi + ... + fij^Zx + '^iZx+i + ... + 'y^Zk+l ( 3 - 2 ) 


• 1 is the multiplicative identity, or real unit. 

Its representation is the unit matrix. 

• /3 is a bireal unit, i. e. it squares to + 1 ; 0 ^ = 1 . 

Its real representation is a symmetric, traceless matrix. 

• 7 is a complex unit, i. e. it squares to —1: 7^ = —1. 

Its real representation is an anti-symmetric (and therefore traceless) matrix. 

The real representation of real, bireal and complex units is orthogonal, i. e. the 
transpose of a representative matrix is its inverse. 

Now let 

7A = IZq -|- fiiZi f3j^ZM — 7i-^M+i... — 'Jn^m+n (3-3) 

the transpose of Z. 


Also, let 


3?Z = Zo 


( 3 . 4 ) 


the real part of Z. 

If Z is a d X d matrix, then 

3 f?Z = iTrZ ( 3 . 5 ) 

since all imaginary (bireal and complex) units have trace 0. 

It is easy to verify that 


(A|B) = 3 fJ(ATB) ( 3 . 6 ) 

is a scalar product. What’s more, all the units form an orthogonal basis with respect to 
this scalar product. Therefore, the remaining coefficients of Z are 


Zk = iR (/31Z) (l<k<K) 


Z*-+I = »(7;z) (1</<L) 


( 3 , 7 ) 



Let us distinguish between an abstract algebra and its matrix representation. An 
element of an abstract algebra also looks like eqn. fl 3 . 2 [) . The addition of two algebra 
elements given in this form is straightforward; multiplication is specihed with a 
multiplication table of the units The matrix representation of the units 

together with matrix multiplication must reproduce this table. 


3.2 Example: the complex numbers 

For example, take the algebra of complex numbers Z = IZq + 72’i, with Zq and Zi real. 
The multiplication table is To save space, let’s truncate this to 



1 7 

1 

7 

1 7 

7 -1 


1 

7 

7 

-1 


The multiplication of two complex numbers is the familiar 

A • B = (lAo + jAi){lBo + jBi) = l{AoBo - AiBi) + jiAoBi + AiBq) 

The real representation of Z is a 2 x 2 matrix, e. g. 

Z = lZo + 7 ^i= ( ( 3 . 8 ) 

Matrix multiplication now reproduces the abstract multiplication law. 

3.3 Clifford algebras 

Here is a list of the simplest associative algebras. 

• The only algebra with 1 unit is the real numbers. 

• There are two algebras with 2 units: the bireal (or split-complex) numbers 
Z = IX -|- I 3 Y, and the complex numbers Z = IX -|- 'yY. 

• There are no algebras with 3 units. 

• Algebras with 4 units are the Cockle (or split-) quaternions 

Z = IZq -|- PiZi -|- ( 32^2 + 7.^3, the tessarines Z = IZi + (3Z2 + + 72 -^ 4 ) 

the Hamilton quaternions: Z = IZq + 'JiZi + ^2^2 + 73-^3. 

How can we construct associative algebras in a systematic way? There are basically two 
ways to do this: by forming the tensor product of two algebras, and with a scheme due 
to W. K. Clifford. I will present the latter, since it automatically classihes the algebra. 
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0.0 Clifford algebras over a M are denoted by The units and multiplication 

table are obtained as follows. 

1 . Start with the real unit 1 and the N = p + q gener¬ 
ator units /3i ... /3p, 7i • • • 7^- The symmetric part of 
the multiplication table of these units is given to the 
right. Any two different generators anti-commute! 

2 . Dehne N{N — l )/2 new units as product of two different generators: 

Iq+l = = -/32/3i, lq+2 = I1I2 = “7271, • • • , f^p+l = f^lll = -7l/3l, ■ ■ • 

3 . Dehne N{N — 1 )(A^ — 2)/6 new units as product of three different generators. 



1 

Pi 

h 

1 

1 

Pi 

In 

Pk 

^k 

dkP 

0 

Ik 

Ik 

0 

^mn 1 


N. Dehne the last new unit u) = I 3 i... / 3 p 7 i... 7^. 

The full multiplication table of all the 2^ units now follows from the symmetric 
multiplication table of the N generator units. 

3.4 Matrix representations 

Here is a table of irreducible matrix representations of Clihord algebras Clp^q with 
p, g < 4 . 



0 

1 

2 

3 

4 

0 

M 

c 

H 

El©]! 

e(2) 

1 

M © M 

M(2) 

C( 2 ) 

H(2) 

e(2) © H(2) 

2 

M(2) 

M(2) ©M(2) 

M( 4 ) 

C( 4 ) 

H( 4 ) 

3 

C( 2 ) 

M( 4 ) 

M( 4 ) © M( 4 ) 

M(8) 

C(8) 

4 

e(2) 

C( 4 ) 

M(8) 

M(8) ©M(8) 

M( 16 ) 


It turns out that every Clihord algebra is isomorphic to 

• a full d X d real matrix algebra or 

• a full d X d complex matrix algebra C{d), or 

• a full d X d quaternionic matrix algebra ]H[((i), or 

• to the direct sum of two such algebras. 

It also appears that some Clihord algebras with diherent p, q, but the same overall 
number of generators p + q, are isomorphic, i. e. their multiplication tables can be 
brought to coincide after re-numbering the units (sign-changes may be necessary, too). 

We are interested in the algebra of the real 2 x 2 matrices M( 2 ), and of the real 4 x 4 
matrices R( 4 ). Both algebras can be constructed with Clihord’s scheme: 

M( 2 ) ^ CkpiR) = M( 4 ) ^ Ch,i{R) = Ck,2{^) 

The multiplication table and representation of these two algebras is given in sections 4 
and 5 . For the other Clihord algebras see ref. [ 5 ]. 


10 















4 One particle degree of freedom: linear optics with real Pauli 
matrices 

4.1 Cockle quaternions 

The real 2x2 matrices M(2) are a representation of the Clifford-Algebras < 7 / 2,0 (I^) and 
C/ 1,1 (M). This is the algebra of the Cockle quaternions (or split-quaternions). 

Units of C/ 2 , 0 (®): f^i, l32^"f = 

Units of C/i 4 (M): 1, (3^,^, /32 = /3i7 


Multiplication table 


1 

/3i 

(32 

7 

/3i 

1 

7 

to 

(^2 

-7 

1 

-/3i 

7 

~ f ^2 

/3i 

-1 


Blue: anti-commuting products. 


General Cockle quaternion and its real representation 


Z — IZq + PiZi -\- (32^2 + 7-^3 — 


Zq + Zi Z2 + Z^ 
Z 2 — ^3 Zq — Z\ 


(4.1) 


Remarks 

• There are 1 real (1), 2 bireal (/ 3 i, /32) and 1 complex (7) units. 
They are represented by the ’’real Pauli matrices” ao, as, cxi, ia2. 

• The general symmetric 2 x 2 matrix is S = ISq -|- / 3 iSi -|- [32^2- 

• The general skew-symmetric 2 x 2 matrix is A = 7A3 . 

• The general diagonal matrix is D = IHq + f 3 iDi . 

• 7 is the symplectic form. 


Determinant and Eigenvalues 


det Z = ZZ = {Zof - Zl - Zl + 


A± = Zo ± ^Jzl + Zl-Zl (4.2) 


Inverse matrix 


Z _ IZo — (3iZi — ( 32 Z 2 — 7^3 
7^ ~ Zl - Zl - Zl + Zl 


( 4 . 3 ) 
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Vector notation 


In the following it will be helpful to combine the part ( 3 iZi + /32-^2 of Z to a vector: 

Z = IZo + P ■ Z + 7Z3 where ^ ^ ^ ^ 

4.2 The beam matrix 

The general beam matrix is 

S = ISo + ^ ■ S - So - El ) 

The beam emittance is 

e = Vdet (S) = ^S 2 -S 2 ( 4 . 6 ) 

A beam matrix in normal form has vanishing S, and is 

S = (4.7) 

4.3 Symplectic transformations 

The symplectic condition R7R''' = 7 translates to det R = (Ro)^ — R\ — R2 + R\ = 1 - 

This formula is not really helpful for constructing symplectic transformations. However, 
it is easy to write down the logarithm of a symplectic transformation: 

In R = r = ^ ■ r + 'yr^, ( 4 . 8 ) 

This is a 7-symmetric matrix, see eqn. (ITBD . 


4.4 Elementary symplectic transformations 

The RHS of eqn. fl 2 . 6 p is a linear combination of the ’’units” (3 ■ e (squares to + 1 , since 
e is a unit vector) and 7 (squares to — 1 ). These units generate the ’’elementary 
transformations”: 

= 1 coshy +/3 • esinhy ( 4 . 9 ) 

= 1 cos V’+ Tsin-^ ( 4 - 10 ) 


Boosts: 

Rotations: 


R-/3 = exp [J 3 ■ ex 

R.y = exp ( 7 -^) 
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Boosting S 


To describe the action of a boost along the unit vector e = (ci, 62 ), let 

S = eS|| + nEj_ 

where 

S|| = e • S Sj_ = n-S = eAS 

The action of a Boost on S is : 

S' = R/ 3 SRT 

= exp (^0 ■ ex) (lT,o + ^ ■ s) exp (^0 ■ ex) 

= 1 (cosh(2x)So + sinh(2x)S||) + 0 ■ [e (cosh(2x)S|| + sinh(2x)So) + 

The components of the transformed 
beam matrix are 

Eg = cosh(2x)So + sinh( 2 x)S|| 

Sj| = cosh(2x)S|| +sinh(2x)So 
S'^ = Ex (4.13) 

Figure 1 : 

R /3 = exp f3 ■ ex boosts a phase vector 
X on a hyperbola with e as the conju¬ 
gate axis, and boosts the vector part of 
a beam matrix S = IEq -|- /3 ■ E in the 
direction of e. 


and E_l 


R = exp (^/3 • exj ^ 

leave the direction of E invariant. If the boost angle is given by 

ta,ih(2x) = (4,16) 

then E' vanishes. 


(4.14) 

(4.15) 


Boost: 5c=0.50; i:q= 1 .50, i:'o=1.47 



The boost invariants are 

E^-E^ 

Boosts parallel to E, i. e. of the form 



(4.11) 

(4.12) 
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Rotating S 


The action of a rotation on S is (see also figure 2): 


S' = R^SRT 

= exp ( 7 ^/^) (lEo + ,9 ■ s) exp 
= ISo + 0 ■ S cos(2-^) + eS sin(2-^) 


Here, e is the skew-symmetric tensor (it looks like 7 ), and eS = (S 2 , —Si)'''. 


Rotation: \|/=0.79; i:q=i:'q=1 .50 



The components of the transformed 
beam matrix are 

S' = So (4.17) 

S' = S cos(2'^) -|- eS sin(2'^) 

Figure 2 : 

R^ = exp 7'0 rotates a phase vector X 
clockwise by the angle 'll), and rotates 
the vector part of a beam matrix S = 
ISq + f3 ■ T, clockwise by the angle 2 'if), 
while leaving the scalar part invariant. 


The rotational invariants are 

So and S^ (4-18) 


4.5 Diagonalizing the beam matrix 

The beam matrix is diagonal when S 12 = S 21 = 0 (standard matrix basis), or when 
S 2 = 0 (real Pauli matrix basis). The two phase space coordinates x and x' are then 
uncorrelated, or decoupled. 

To diagonalize a beam matrix S, decompose S = ISo -|- /3 ■ S. Rotate S with 

S2 

R^ = exp (j'lpi) tan( 2 ?/>i) = — 

The new beam matrix is diagonal: S' = R.^SRt = ISq + /^i^n 

The diagonal form is preserved under the group of boosts R^ = exp {(3ix) (free param¬ 
eter x). 
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4.6 Normalizing the beam matrix 


The beam matrix is normal when S 12 = S 21 = 0 and Sn = S 22 = e (standard matrix 
basis), or when Si = S 2 = 0 (real Pauli matrix basis). 

To normalize a beam matrix S, hrst diagonalize it to S' = ISg + /3iS'j^, see section 4.5. 
Then boost S' with 

S' 

R/3 = exp (/3ix) tanh(2x) = 

The new beam matrix is in normal form: S" = R^S'R^ = le. 

The normal form is preserved under the group of rotations Rg, = exp ( 7 "^) 

(free parameter x)- 

The normalization can also be done without the intermediate diagonalization, as follows. 
To normalize a beam matrix S, decompose S = ISq + /3 ■ S. Boost S with 

R /3 = exp ($-ex) A tanh( 2 x) = -|p- 

^ '' |S| ^0 

The new beam matrix is in normal form: S' = R^SR^ = le. 


4.7 Invariance group of the beam matrix 

Which transformations I leave the beam matrix S invariant: ISIt = S ? 

Let N be a normalizing transformation to a beam matrix S, i. e. S = NleN'''; the 
beam matrix is invariant under the group of transformations I {ip) = N exp {'yip) 

(free parameter ip). 
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5 Two particle degrees of freedom: linear optics with real Dirac 
matrices 


5.1 The Clifford algebra 


= 


The real 4x4 matrices M( 2 ) are a representation of the Clifford-Algebra Cl^^i 
Cl 2 , 2 (M), see section 3.4. The 16 units of this algebra can be grouped into 4 ’’scalars”: 

1 , 7^ and 4 ”3-vectors”: Ck ^ C Pi ^ Pi -t Pi with k = 1,2,3. 


Units of 

1 Pln^ Pk= PWXk = hklmPlPln Pk = CkJ^ l‘^ = -Plp2pl 7 ^ = 7 ^ 7 ^ 
Units of C/ 2,2 (IK): 

1 Pi = Ph\pl = Ph^Pl = -Pli\Pl = -Ph\Ci = pIpW = 7V 

Pi = PlPh^Pl = PIpW, C 2 = PWi\C^ = PWi^ P\ = -PlPWi^ 


Multiplication table 


1 


Cl 

P? 


^mno'y 

PT 

^mnCl ^mnof^l 

Ck 

PI 

—^kP — ^kliCi 

—^klln ~ ^kliP7 

Pk 

^mnof^k 

~^kllm ~~ ^kliPT 

^kl^mn^ ^kl^mno^o 

• • “1“ ^kli^mnCi ^kli^mnof^i 


Blue: anti-commuting products. 


General C/34-sedenion and its representation 


Sedenion means ’’number with 16 components”. These components are the coefficients 
that come with the units. Like the units, the coefficients can be grouped into 4 
’’scalars”: Z^, and 4 ”3-vectors”: Z° ^ Z°, Z^ ^ Z^ Z| ^ Z^, Z| ^ Z^, 

with A; = 1, 2, 3. In products of units and/or coefficients with indices, summation over 
equal indices is required. 


Z = 1 Z° + CkZl + + PkZk = 1^0 + 6 + 7 ™^o" + P ""-= ( 5 - 1 ) 


/ Z 0 + Z/ + Z| + Z| 
Z? - Zl - Z| + Z| 

7O I 72 73 

Z/o -|- Zjo — Z/n — Z/i 

V zl-zl + zl-zl 


zl + zl - zl + zl 
zl + zl - zl - z| 
zl + zl + zf - zl 
zl + zl + Zq + zf 


zl + Z3I + zl - zl 

Zg -|- zl -\- z! -\- Zg 

70 71 I 72 73 

Z/n — \ -^9 — 

z?-zi + z| + z| 


-Z 3 O - zl + zl + zl \ 
-zl + Z 3 I - zl + zl 
Z? + Zgl + Z| + Z| 

Z0_^l_^2 + ^3 


See also appendix B.l. 
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Remarks 


• There are 1 real (1), 9 bireal (/3™) and 6 complex 7 ™) units {k,m = 1,2,3). 
They are represented by ’’real Dirac matrices”. 

• The general symmetric 4x4 matrix is S = ISq + . 

• The general skew-symmetric 4x4 matrix is A = Ck^k + . 

• The general diagonal matrix is D = IDq -|- f3\Dl + 0lDl + fSlD^ . 

• 7 ^ is the symplectic form. 

• The 16 components of the number Z are arranged as a 4 x 4 matrix, not as a 
vector. Don’t confuse this matrix with the representative matrix (Zf^^)! 

Component matrix of Z; 

/zg zg zg zg\ 

yO yl y’l yZ 

(y\\ _ ^1 

) ^0 ^2 ^3 

y zl z] z\ z\ y 


Representative matrix of Z: 


(Z,.) = 


/ ^11 

Zl2 

Zl3 

Zi 4 

\ 

Z21 

Z22 

Z23 

■^24 


Z31 

Z32 

Z33 

■^34 


\ Zil 

Zi2 

Zi3 

Z44 

/ 


Determinants 


The anti-symmetric 4x4 matrix A = -l- 7 ™'A™ has the determinant 

det A = [{Air+{Air+{Air - {Air - i^ir - (a^)^]' 

Let A = S 7 iS = CfcA 0 + 7 -A- 

where Aj = (S°) V SL Ag = 2SL Al = 2 E^-E^ 

and A° = 2(^S2 aS3-S°Si) 

The symmetric 4x4 matrix S — ISg + /3’"- S’" has the determinant 

det S = Vdet (S7'E) = VdetA = - ri"- ri" 


(5.2) 


(5.3) 


Let A = Z 7 ‘Zt = + 7 ’"Ag‘ 

where AJ = ZgZg + ZgZg - ZgZg - ZjZg 

K = 2 (ZgZg - ZgZg) Ag = 2 (ZgZg A ZgZg) 
and A“ = 2 izgZ” - ZgZ‘ + Z“a Z' + ZgZ^ - ZgZ® + ZV zA 

The general 4x4 matrix IZ^ ^ -|- 7 ™'Z™ -|- 0"^- Z^ has the determinant 

det Z = v^det (Z 7 IZT) = Vdet A = (Aj)^ + {Air + (Aj))^ - A°- A° (5.4) 


(Choosing 7 ^ instead of 7 ^, 7 ^, Cit C2 C3 for constructing a skew-symmetric matrix A 
out of S or Z is arbitrary. But only A = S 7 ^S transforms in the same way as S, to 

wit A' = RART.) 
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Inverse matrix 


The anti-symmetric 4x4 matrix A = + 7 ™"A™ has the inverse 

A-i_ + 

~ {Alf + (A2)2 + (yl3)2 - (yl?)2 - - (AO)^ 

The symmetric 4x4 matrix S = ISg -|- /3”^- has the inverse 

= 7 ^SA~^ where A = S 7 ^S 

For the components of A in terms of S see the paragraph Determinants. 

The general 4x4 matrix IZq + ^ Z™ has the inverse 

where A = Z 7 ^Zt 

For the components of A in terms of Z see the paragraph Determinants. 


(5.5) 


(5.6) 


(5.7) 


5.2 The beam matrix 

The general beam matrix is 

/ s[] + sj + Si + S3 -si + S3 -S3 + si -si + s? \ 

-si + si sg + s{-si-si si + sf sf + si 

-S3 + Si si + s? Si-S} + Si-S3 Si + S3 

I vl _i_ v2 ^*3 _i_ vl ^^2 _i_ ^*3 vO vl v2 i v3 , 

\ ~^2 "T “T ^3 ^3 -h Zj2 Zjq — Zj;^ — Zj 2 -h Zj3 / 

The two emittances of S appear in the characteristic polynomial of S 7 ^, see section 2.3: 

det (S 7 ^ - Al) = (A^ + ei) (A^ + ejj) (5.9) 

On the other hand we have 


det (S 7 ^ - Al) = A/det (A -h = (AJ -h A^)^ -h (A^)^ -h (Ag)^ - A^- A^ (5.10) 
with A = S 7 ^S given in the paragraph Determinants. The comparison gives 

r,J, = K±^JA«■A«-(Alf-^A:l;f (5-11) 

Like the emittances, the quantities Aj and A^- — (Ai)^ — (A)])^ are transformation 

invariants. 
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5.3 3D visualization and numerical example 

Here is a stereo view of a random beam coefficient matrix: yon shonld look at the left 
fignre with your left eye, and at the right hgure with your right eye. The beam matrix 
consists of a scalar Sq, shown as a gray bar along the upright 1-axis, and three vectors 
shown in primary colors. The scalar products of the three beam vectors are 
given in the figure heading: (R,G) means ’’red vector dot green vector” etc. The several 
decoupling transformations in the later sections will all be applied to this initial beam. 
The visualization is often helpful for hnding a decoupling procedure. 


(R,G) = -1.43, (R.B) = -4.41, {G,B) = 2.18 



(R,G) = -1.43, (R,B) = -4.41, (G,B) = 2.18 



Figure 3: 


Random initial 

beam 

with coefficients 

S': 


/ 4.4 

0 0 

0 \ 


0 

1.5 1.6 

0.9 


0 

1.5 0.7 

0.9 


V 0 

1.7 1.6 

2.6 y 



1 -axis upright 
.3 2 -axis to your left 
3-axis to your right 


5.4 Symplectic transformations 

The symplectic condition A = = 7 translates to 

AJ = RlRl + RlRl - RIr:^ - R^Rl = det R = 1 

Ai = 2 (RlRl - RjRy =0 Al = 2 (KlRl + RJr;) = 0 

A“ = 2 + R“a R' + R^S‘ - «=>) = 0 

It can be shown that the conditions Ag = Ag = 0 split up into 

pO p2 _ pO p3 _ pi p2 _ pi p3 _ n 

and that the conditions A° = 0 splits up into 

+ iP A R^ = RlR^ - R^R^ + R° A R^ = 0 

These formulas are not very helpful for constructing symplectic transformations. How¬ 
ever, it is easy to write down the logarithm of a symplectic transformation: 

In R = r = ^ ■ r ° 7 ^rQ -|- ^9^ • r ^ • r ^ (5.12) 

This is a 7 -symmetric matrix, see eqn. 02.611 . 
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5.5 Elementary symplectic transformations 


The RHS of eqn. fl5.12p is a linear combination of the four units C ' (squares to — 1), 
7 (squares to —1), f3‘^- e'^ (squares to +1) /3^ - (squares to +1) , where is a unit 
vector. These units generate the ’’elementary transformations”: 


C-Rot at ions: 

Rc 

= exp 

(c • #) 

= 1 cos "0 + C ■ e sin 

(5.13) 

7 -Rotations: 

Rg, 

= exp 

(tV) 

= 1 COS0 -f 7 ^ sin0 

(5.14) 

/3^-Boosts: 

R/32 

= exp 

ex) 

= 1 cosh X + • e sinh x 

(5.15) 

/3^-Boosts: 

R-/33 

= exp 

ex) 

= 1 cosh X + /3^ • e sinh x 

(5.16) 


Below I give the transformed beam matrix S' = RSRt under the above elementary 
transformations. The calculation is laborious, details are given in Appendix B.2. 

Earlier we have found that the invariants of S under the 10 parameter group of 
symplectic transformations are the two emittances e/, e//. Under the 1-parameter 
subgroup of elementary elementary transformations, S has many additional invariants. 
These will be listed together with the transformation laws. 


(^-Rotation 


The action of a (^-rotation on S is: 


R^ERJ 







exp 

(«■ 

cip 

)(eS + 

0^- 

S™) exp 

'-C • eiA 

^^0 + 









+ 

COS (2yi) 

Si 

— sin (2'0) 

e AE^ 

+ 



+ 

cos (2'0) 

Ei 

— sin {2il>) 

e AE^ 

+ 



+ 

cos (2-^) 

si 

— sin (2'0) 

e AE^ 



(5.17) 




Figure 4: 

A ^-transformation 
leaves Eg unchanged 
and rotates S^, T? and 
cw around e by 2'0. 

(e = (l,l,l)V\/3) 
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The transformation invariants are 


Sq all scalar prodncts S*-E™ where Z,m = 1,2,3 Em Em E| (5.18) 


7 -Rotation 


The action of a 7 ^-rotation on S is: 


S' 


exp (7V) (lS° + r- S'") exp (-7V) 
lE° + ^Sgi + 




cos ( 20 ) E^ -|- sin ( 20 ) E^ 


+ 


— sin (20) -I- cos (20) 


(5.19) 



Fignre 5: 

A 7 -transformation 
leaves Eq and E^ nn- 
changed and rotates 
E^ & along the in¬ 
variant ellipse spanned 
by these two vectors. 


The transformation invariants are 

E° S^AS^ (S2)2 + (S3)2 (S^. S^)2 + (SS S^)2 (5.20) 

To rotate E'^ and E'^ into main axes of the invariant ellipse (E'^- E'^ = 0), choose the 
rotation angle according to 


tan(40) 


2S2. S3 
(E2)2 - (S3)2 


(5.21) 
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/3^-Boost 


The action of a /3^-boost on S is: 

= Rfl 2 SR ,^2 


■132 

= exp ex) + 0"^- E”^) exp ex 

= 1 [cosh (2x) E[^ + sinh (2x) E|] + 

^11 + sinh(2x)e A E^ + cosh(2x)E[L + 
e (sinh (2x) Eq + cosh (2x) Ey) + + 

■ eL^ + cosh (2x) E[]_ — sinh (2x) e A 


(5.22) 


where Ei|^ = e • E™ and E™ = E™' — eEi, 



Figure 6: 

A /^^-transformation 
mixes Eq and E|, and 

mixes E)|_ and E[[. 
(e-'=(l,l,l)V\/3) 


The transformation invariants are 

Em 


(E^)^ - (E^)^ (E^)^ - (E^)^ 


S^- 

e^ae^-e^e^V- ("eAe^"' 


A /3^-boost also allows to supress E^. This is achieved by setting 


e = 


E 2 

|E2| 


tanh(2x) = 




A /3^-boost allows to adjust S2- while keeping S^- constant. 
To make S'2- S'^ = 0, set 

S 2 .S 1 


e = 


E2aE3-E°E1 


tanh(2x) = 


|E2a E3 - E[]E1| |E2AE3-E[jEi| 

A /3^-boost allows to adjust E2- E^ while keeping E^- E^ constant. 

To make S'2 - S'^ = 0, set 


e = 


S^AS2-SgS3 

|SiAS2-E[]S3| 


tanh(2x) = 


S2.E3 


|SiAS2-E[jE3| 


(5.23) 


(5.24) 


(5.25) 


(5.26) 
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/3^-Boost 


The action of a /3^-boost on S is: 


'/33' 

= exp (^^3. exp 

= 1 [cosh (2x) So + sinh (2x) S|] + 


• ex 


Sill — sinh(2x)e A S^ + cosh(2x)S[|_ + 
Sly + cosh (2x) S^ + sinh (2x) e A S^ + 
S^- [e (sitih (2x) ES + cosh (2x) Sf) + 2= 


(5.27) 




Figure 7: 

A /^^-transformation 
mixes Sq and S|, and 

mixes S^ and S^. 

(e = (l,l,l)V\/3) 


The transformation invariants are: 

S^ Sn 


Si 


SA 


(S°)2 - (E^)^ (S^)^ - (S^) 


^352 


^ 2\2 


s"as"-s"s' 


SS 


A /3^-boost also allows to suppress S^. This is achieved by setting 

^ S3 


e = 


|S3| 


tanh(2x) = 


A /3^-boost allows to adjust S^- E^ while keeping E^- E^ constant. 
To make E'^- E'^ = 0, set 


S^AS^-EgSi 

|S2aS3-E[]S1 


tanh(2x) 


S3-S1 

|S2AE3-E[jEi 


A /3^-boost allows to adjust S3- while keeping S^- constant. 
To make S'3- S'^ = 0, set 


E3a E^ - E[]E2 
|E3aE1 -E[]E2| 


tanh(2x) 


E3-E2 

|E3aE 1 -E[jE2| 


(5.28) 


(5.29) 


(5.30) 


(5.31) 
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5.6 Decoupling two pairs of phase space coordinates 

Decoupling two complementary sets A and A of phase space coordinates means, making 
the correlation between A and A vanish. We take the 4 phase space coordinates to be 
{Xi, X 2 , X 3 , X 4 } = {x, x', y, y'} {=AuA). 


Decoupling x,x' from y,y' 


This means bringing the representative matrix and component matrix into the form 


^ Sii S 12 0 0 ^ 

S 21 S 22 0 0 

0 0 S 33 S 34 

\ 0 0 S 43 S 44 y 


/ s[j 0 0 0 \ 

0 sj 0 0 

0 0 S3 

V 0 0 S3 S3 y 


To decouple x,x' from y,y’, proceed as follows (steps 1 and 2 can be interchanged): 

1. Decompose S = ISg + S™'. 

Boost S^- S^ to 0 with S' = R^ 2 SR^ 2 ) where 


R/32 = exp /3 - ex 


e = 


S2AS3-S[jS3 

|S2aS3-S0S1| 


tanh( 2 x) = 


ST Si 


|S2aS3-S[]S1| 


2. Rename S' —)■ S and decompose S = ISq + S'". 

Boost S3- Si to 0 (keeping S^- Si = 0) with S' = R^ 3 SR^ 3 , where 


R /33 = exp /3 - ex 


e = 


S2AS3-S[jSi 

|S2aS3-S0S1| 


tanh( 2 x) = 


S3. Si 


|S2aS3-S[]S1| 


3. Rename S' ^ S and decompose S = ISq + S'". 

Rotate si parallel to the 1-axis with S' = R^SRJ, where 


R^ = exp 


e = 


( 0 ,si-si)T 

V(E>) 2 +(E‘P 


tan( 2 y’) = 


y(E|)^ + (Ej)^ 

s! 


The decoupling is preserved under the transformation group R = exp [Ci^P + 7^0 
+f3lxi + 01 X 2 + + / 33 X 4 ] (free parameters 0, 0, Xi, X 2 , Xs, X4). 
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{R,G) = 0.00, (R,B) = -4.41, {G,B) = 0.57 (R,G) = 0.00, (R.B) = -4.41, (G,B) = 0.57 




Figure 8: 


Dec. X 

x' 

y,y' 

step 

S3 

ffi s 

2 



Coefficient 

matrix S^ 

/ 

4.3 

0 

0 



0 

1.3 

It—1 

0.7 


0 

1.6 

0.6 

0.8 

V 

0 

1.8 

0.9 

2.7 / 


(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.46 



(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.46 



Figure 9: 

Dec. x, x' I y, y', step 2: 
E2, S3 

Coefficient matrix S^ = 



/ 3.5 

0 

0 



0 

I.l 

1.6 

0.5 


0 

0.3 

0.5 

1.2 

3 

^ 0 

2.1 

0.7 

0.4 / 


(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.46 (R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.46 



Figure 10: 

Dec. X, x' I y, y', step 3: 
Siffi S2, S^ffi S3 
S^ II 1-axis 


Coefficient matrix S^ = 


/ 3.5 0 0 0 \ 

0 2.3 0 0 

0 0 0.6 I.l 

\ 0 0 1.7 0.7 
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Decoupling x, y from x', y' 


This means bringing the representative matrix and component matrix into the form 


S = 


/ Ell 0 Si3 0 \ 

0 S 22 0 Si4 

S34 0 E33 0 

V 0 S41 0 S44 / 


(s^) = 


/ Eg 0 0 0 \ 

0 Ej 0 Ef 

0 0 E2 0 

Vo Eg 0 Eg y 


To decouple x, y from x', y', proceed as follows: 

1. Decompose S = lEg + 0 ^- S”^. 

Boost E^- E^ to 0 with S' = R^ 2 SR ^25 where 


R 


./32 = exp 


(/3= 


• ex 


e = 


E2aE3-E[]E^ 

|E2aE3-E[]E1| 


tanh(2x) = 


ET El 


|E2aE3-E[]E1| 


2. Rename S' —)■ S and decompose S = IEq + ^9™'- E™'. 

Boost E3- E^ to 0 (keeping E^- E^ = 0) with S' = R^ 3 SR^g, where 


R «3 = exp 


(/3= 


• ex 


e = 


S3aE1-E°E2 
|S3aE1 -S[]S2| 


tanh(2x) = 


S3-E2 


|E3aE1 -E[]E2| 


3. Rename S' —)■ S and decompose S = IEq + / 3 ^- E™'. 
Rotate E^ parallel to the 2-axis with S' = R^SRJ, where 


R^ = exp 


e = 


(-sj.o.sjr 

V(E?)2 + (Ei)2 


tan(2'0) = 


V(e;)^ + (Ej)^ 
Si 


The decoupling is preserved under the transformation group R = exp [<^ 2 '*/' + 01x1 
+0iX2 + 0lx^\ (free parameters •0, Xi, X 2 , Xs)- 


{R,G) = 0.00, (R,B) = -4.41, (G,B) = 0.57 (R,G) = 0.00, (R.B) = -4.41, (G,B) = 0.57 



Figure 11: 

Dec. x,x' I y, ?/', step 1: 
E^T E2 

Coefficient matrix E^ = 


/ 4.3 

0 

0 


0 

1.3 

1.9 

0.7 

0 

1.6 

0.6 

0.8 

V 0 

1.8 

0.9 

2.7 / 
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{R,G) = 0.00, (R,B) = -4.34, {G,B) = 0.00 


(R,G) = 0.00, (R,B) = -4.34, (G,B) = 0.00 



Figure 12: 

Dec. X, x' I y, y', step 2: 
S3 


Coefficient matrix S^ = 



/ 

4.3 

0 

0 




0 

1.3 

1.8 

1.3 



0 

1.5 

0.2 

0.9 

3 

V 

0 

1.7 

I.l 

2.3 / 


{R,G) = 0.00, (R,B) = -4.34, (G,B) = 0.00 (R,G) = 0.00, (R,B) = -4.34, (G,B) = 0.00 



Figure 13: 

Dec. X, x' I y, y', step 3: 
S^ffi S2, S^ffi S3 
S^ II 2-axis 

Coefficient matrix S^ = 


/ 

4.3 

0 

0 

0 \ 


0 

2.5 

0 

0.6 


0 

0 

2.2 

0 

V 

0 

1.0 

0 

2.7 / 
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Decoupling x,y' from x',y 


This means bringing the representative matrix and component matrix into the form 


/ Sn 

0 

0 

Ei4 

\ 

0 

^22 

^23 

0 


0 

^32 

^33 

0 


V ^41 

0 

0 

E 44 

/ 


/ EO 0 0 0 \ 

0 Ej S3 0 

0 S2 S2 0 

\ 0 0 0 E| / 


To decouple x, y' from x', y, proceed as follows: 

1. Decompose S = ISg + 0’^- S'". 

Boost E^- E^ to 0 with S' = R^aSR^g, where 


R 


■/33 = exp 


(/3= 


• ex 


e = 


E2aE3-E[]E^ 

|E2aE3-E[]E1| 


tanh(2x) = 


ED S3 


|E2aE3-E[]E1| 


2. Rename S' —)■ S and decompose S = ISq + S'". 

Boost E^- E3 to 0 (keeping E^- E3 = 0) with S' = R^ 2 SR^ 2 ; where 


R ^2 = exp [p" 


■ ex 


e = 


E^AE2-E°E3 

|S1aE2-S[]S3| 


tanh(2x) = 


ets3 


|E1aS2-E[]E3| 


3. Rename S' —)■ S and decompose S = IEq + /3'"- S'". 
Rotate E3 parallel to the 3-axis with S' = R^SRJ, where 


R^ = exp 


e = 


(s|.-s;,o)T 

/(sfFRElp 


tan(2'0) = 


+ (Sj)^ 

Si 


The decoupling is preserved under the transformation group R = exp [<^ 2 '*/' + 01x1 
+0iX2 + 0lx^\ (free parameters •0, Xi, X 2 , Xs)- 


{R,G) = -1.43, {R,B) = 0.00, (G,B) = 0.26 (R,G) = -1.43, (R,B) = 0.00, (G,B) = 0.26 



Figure 14: 

Dec. x, x' I y, 2 /', step 1: 
E^T S3 

Coefficient matrix S^ = 


/ 3.6 

0 

0 

0 \ 

0 

1.3 

1.2 

0.3 

0 

0.2 

0.4 

1.3 

V 0 

2.0 

1.5 

0.3 / 
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{R,G) = -1.26, {R,B) = 0.00, {G,B) = 0.00 


(R,G) = -1.26, (R,B) = 0.00, (G,B) = 0.00 




Figure 15: 

Dec. X, x' I y, y', step 2: 

E3, S3 

Coefficient matrix S^ = 



/ 3.2 

0 

0 

0 \ 


0 

1.2 

0.1 

0.6 


0 

0.2 

0.9 

0.2 

3 

V 0 

1.7 

0.5 

0.4; 


{R,G) = -1.26, {R,B) = 0.00, (G,B) = 0.00 (R,G) = -1.26, (R,B) = 0.00, (G,B) = 0.00 



Figure 16: 

Dec. X, x' I y, y', step 3: 
S^ffi S3, S^ffi S3 
S^ II 1-axis 


Coefficient matrix S^ = 
/ 3.2 0 0 0 \ 

0 2.1 0.4 0 

0 0.5 1.0 0 

\ 0 0 0 0.7 / 


5.7 Diagonalizing the beam matrix 


Diagonalizing the beam matrix means mutually decoupling all 4 phase space 
coordinates. The beam component matrix is then also diagonal, i. e. it has S^, S^ and 
S3 aligned to the 1-, 2- and 3-axis respectively. 


/ Sn 0 0 0 \ 

0 S22 0 0 

0 0 S33 0 

y 0 0 0 S44 J 




/ Sg 0 0 0 \ 

0 Sj 0 0 

0 0 si 0 

Vo 0 0 S3; 
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To diagonalize a beam matrix, first transform the beam matrix to block-diagonal form, 
see steps 1 to 3 in section 5.6. Then continue as follows. 

4. Rename S' ^ S and decompose S = ISq -|- 0’^- S'”. 

Rotate S^- S^ to 0 (keeping S^- S^ = S^- S^ = 0) with S' = R.ySRT —>■ S, where 

2 ST S3 
(S 2)2 - (S3)2 

5. Rename S' ^ S and decompose S = ISq -|- 0’^- S'”. 

The three vectors S^, S^, S^ are now mutually perpendicular. 

Align S^ and S^ to the 2- and 3-axis with S' = R^SRJ —)■ S, where 

T? S3 

Rc = exp (CiV’) tan(2?/^) = ^ = 

The diagonal form is preserved under the transformation group R = exp [(3lx2 + f^lxs] 
(free parameters X 2 ,X 3 )- 


Rq, = exp ( 7 ^ 0 ) tan(40) 


(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 (R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 



Figure 17: 

Diagonalization, step 4: 
S^T S2, S^T S3, 

S^T S3 


Coefficient 

matrix S^ 

/ 3.5 

0 

0 


0 

2.4 

0 

0 

0 

0 

0.3 

1.3 

V 0 

0 

1.8 

0.2 / 


(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 



(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.1 



Figure 18: 

Diagonalization, step 5: 
S^ II 1-axis, S^ II 2-axis, 
S3 II 3-axis 

Coefficient matrix S^ = 
/ 3.5 0 0 0 \ 

0 2.4 0 0 

0 0 1.8 0 

\ 0 0 0 1.3 / 
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Here is an alternative diagonalization recipe that doesn’t go throngh 
block-diagonalization. 


To diagonalize a beam matrix, proceed as follows. 

1. Decompose S = ISg + E™'. 

Boost E^- to 0 with S' = R/ 32 SR^ 2 ) where 


R/32 = exp /3 - ex 


e = 


S^AS^-EgS^ 

|S2aS3-E[]S1| 


tanh(2x) = 


E 2 . El 


|E2aE3-E[]E1I 


2. Rename S' —)■ S and decompose S = IEq + S'". 

Boost E^- El to 0 (keeping S^- S^ = 0) with S' = R^gSR^g, where 


R/33 = exp /3 - ex 


e = 


S^AS^-SOSi 

|S2aS3-S[]S1| 


tanh(2x) = 


EE El 


|S2aS3-S[]S1| 


3. Rename S' —)■ S and decompose S = ISq + /3'"- S'". 

Rotate S2- S^ to 0 (keeping S2- S^ = S^- S^ = 0) with S' = R.ySRT, where 


Rg, = exp ( 7 V) tan(4(;/)) = 


2 S 2 . S3 


(S2)2 _ (^3)2 


4. Rename S' —)■ S and decompose S = ISg + 0^- S'". The three vectors S^, S2, 
E3 are now perpendicular to each other. Calculate 


Ri = TEVISI' 


= iSVlS^ 


f j f j ' j <‘2i 


To make the last rotation angle as small as possible, choose the two (independent) 
signs such that Tl + T 2 + T| becomes maximally positive. Calculate 


cos(2^) = (7 + r| + r| -1) /2 

n - n 


ei = 


T! - Ti 


62 = 


63 = 


tjj in [0,7r/2] 

rp2 rp\ 

-'l -‘-2 


2sin(2^/>) 2sin(2^/>) 2sin(2^/>) 

Rotate Si, S2, S3 parallel to the 1, 2, 3-axis with S' = R^SRJ, where 


R^ = exp ( C • 6^ 


Comment; The formulas in step 4 for the ^-rotation parameters e and 2%l) were 
obtained as follows. There are 4 possible ways to align S^ to (1, 0, 0)''', S2 to (0,1, 0)''' 
and S3 to (0, 0,1)''' (two signs are free). The idea is to choose the alternative which 
needs the smallest rotation angle. The 3 unit vectors T^, T'^ and that are rotated 
into (1, 0, 0)''', (0,1, 0)''' and (0, 0,1)''' and e are dehned accordingly. The matrix 
T = [T1,T2,T3] performs the inverse of this rotation in M3_ Eulers rotation formula 
expresses T in terms of e and tjj: Tki = Ski cos(2-^) — eieiki sin(2-^) -|- ekei{l — cos(2-^)). 
Follows now Tkk = 1 -l- 2 cos(2-^) and eikiTki = —2ei sin(2-^), and the formulas for e h ip. 
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{R,G) = 0.00, (R,B) = -4.41, {G,B) = 0.57 


(R,G) = 0.00, (R,B) = -4.41, (G,B) = 0.57 




Figure 19: 

Diagonalization, step 1: 


Coefficient matrix 


/ 4.3 

0 

0 

0 \ 

0 

1.3 

1.9 

0.7 

0 

1.6 

0.6 

0.8 

V 0 

1.8 

0.9 

2.7 / 


(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.46 



(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.46 



Figure 20: 

Diagonalization, step 2: 
E\ fiffi S3 

Coefficient matrix S^ = 


/ 3.5 

0 

0 


0 

I.l 

1.6 

0.5 

0 

0.3 

0.5 

1.2 

V 0 

2.1 

0.7 

0.4 / 


(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 



(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 



Figure 21: 

Diagonalization, step 3: 
S^ffi S2, S^ffi E3 
E^ffi E3 

Coefficient matrix S^ = 


/ 3.5 

0 

0 

0 \ 

0 

I.l 

1.6 

0.0 

0 

0.3 

0.1 

1.3 

V 0 

2.1 

0.8 

0.2 / 


(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 



(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 



Figure 22: 

Diagonalization, step 4: 
E^ II 1-axis, S^ II 2-axis 
S3 II 3-axis 


Coefficient 

matrix S^ 

/ 3.5 

0 

0 


0 

2.4 

0 

0 

0 

0 

1.8 

0 

V 0 

0 

0 

1.3 ) 
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5.8 Normalizing the beam matrix 


A beam matrix in normal form has aligned to the 1-axis, and vanishing 


/ ei 0 0 0 \ 

0 e/ 0 0 

0 0 e// 0 

V 0 0 0 ejj J 




/ EO 0 0 0 \ 

0 Ej 0 0 

0 0 0 0 

\ 0 0 0 0 / 


To normalize a beam matrix, hrst transform it to diagonal form, see sections 5.6-7. 
Then rename S' —)■ S and decompose S = lEg -|- f3\Y\ + 

Annihilate E^ and E| with S = where R^c is the scale transformation 



The normal form is preserved under the transformation group R = exp (CiV' + 7^0) 
(free parameters "0, </). 


Comment: When the beam matrix S is diagonal, a symplectic rescaling of the phase 
space axes achieves normalization. This rescaling could be expressed in terms of a 
double-boost: R^c = exp {(3\xi + 01 x 2 ), but the formulas for the two boost angles are 
unnecessarily complicated. It is simpler to construct R^c with the elements of the 
(diagonal) representative matrix. 


(R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 (R,G) = 0.00, (R,B) = 0.00, (G,B) = 0.00 



Figure 23: 
Normalization: 

E^ II 1, E2 = E3 = 0 

Coefficient matrix E^ = 



0 

0 


0 

2 

0 

0 

0 

0 

0 

0 

V 0 

0 

0 

0 / 


5.9 Invariance group of the beam matrix 

Which transformations I leave the beam matrix S invariant: ISIt = S ? 

Let N be a normalizing transformation to a beam matrix S, i. e. 

S = N (IEq -|- /3iEJ) N'T; the beam matrix is invariant under the group of transforma¬ 
tions I (-0, 0) = N exp (CiV' + 00) (free parameters ip, 0. 
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5.10 Alternative transformation formulas 

I haven’t addressed yet the problem of decoupling one phase space coordinate from the 
other three. The reason is, the visual approach used until now didn’t suggest to me a 
way to achieve this kind of decoupling. Instead I had to resort to the following 
’’alternative rotation formulas”. (The ’’alternative boost formulas” are given for 
completeness; they will not be used.) 

The linear transformation formulas for S given in section 5.5 can be written as a 10 x 10 
matrix acting on a 10-component vector or more generally as a 16 x 16 matrix 
acting on a 16-component vector Z^. Matrix formulas of this size are not very alluring. 
However, there is a more compact way to express the transformations: by 4 x 4 matrices 
acting from the left or right on the component matrix Z^. Here is how it is done. 

First we express the action of multiplying Z from left or right with a unit as a matrix 
multiplication of a unit with the coefficient matrix {Z^)\ 

Z' = CfcZ ^ {Z^J = UZ^,) and Z'= ZC, ^ {Z^J = 

Z' = VZ ^ {Z^J =and Z'= ^ {Z^J = {Z^^W 

Z' = (3iZ => {Z^J =-C\Z^,)C^ and Z'= Z(3i => {Z^J = 

From these formulas, the following transformation formulas can be derived. 

(^-Rotation 

Z' = exp (^C • #) Zexp (^-C • ^ {Z^)' = exp ((C^ + 7^)efc^) {Z^)- 

7 -Rotation 

Z'= exp ( 7 V) Zexp (-7V) ^ (Z^)'= exp (-(Ci + 7^)0)- 

Note that Cfc + 7^ = ‘^Jki where 

00 oo\ /ooo 0\ 

00 00 000 -1 
0001 ^“0000 

00-10/ VOlO 0/ 

/3^-Boost 

Z' = exp j z exp j ^ [Z^]' = \[{Z^) - f3[ekei{Z^)0^^ + ... 

• • • I [{Z^) + cosh2x + I [-Cfcefc(Z^)C 2 - 7''efc(^^)7^] sinh2x 

/3^-Boost 

Z' = exp j Z exp j ^ [Z^)' = I [{Z^) - f3{ekei{Z^)0^^ + ... 

• • • I [{Z^) + cosh2x + 4 [-Cfcefc(Z^)C 3 - 7''efc(^^)7^] sinh2x 
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5.11 Decoupling one phase space coordinate 

In the following, s designates a symmetric element of the component matrix, 

and a designates a skew-symmetric element E^ = —E^. 


Decoupling x from x',y,y' 


This means bringing the representative matrix and component matrix into the form 


S = 


( Sn 

0 

0 

0 

\ 

0 

S22 

S23 

S24 


0 

S32 

S33 

S34 


1 0 

S42 

S43 

S44 

/ 


/ 


(S^) = 


0 

0 Ej 


0 

s 


0 

V 0 


0 \ 

s 
s 

n ) 


To decouple x from x', y, y', calculate the 3D polar decomposition of the ’vector-vector- 


part’ of the component matrix: (Ej,) = 0 {T,\) 
El E? E? 


El Ei 

Si 


Si 

Si 


= exp 


2i) 


0 

ea 

-62 


-63 

0 

6l 



Extract 'ifj and e from O; 

cos{2i/j) = (Ol + Ol + Ol — /2 

O? 


6l = 


Oi 


Di 


2 sin( 2 ^/>) 


62 = 




2 sin( 2 ^/>) 


63 = 


ijj in [ 0 , 7 r/ 2 ] 

01 - Qj 

2 sin( 2 ^/>) 


Then S' = exp ■ ei/j) T, exp ( —^ • e^jJ) has x decoupled from x', y, y'. 


The decoupling is preserved under the transformation group R = exp -|- 7 ^) 

+ I (/^3 + /^i) Xi + 01x2 + /^ixs] (free parameters 0 ,Xi,X 2 ,X 3 )- 


(R,G) = -1.43, (R.B) = -4.41, (G,B) = 2.18 



(R,G) = -1.43, (R,B) = -4.41, (G,B) = 2.18 



Figure 24: 

Beam with x decoupled. 
Coefficient matrix E^ = 


/ 4.4 

0 

0 


0 

2.6 

0.2 

0.8 

0 

0.2 

2.3 

0.4 

V 0 

0.8 

0.4 

2.8 / 
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Decoupling x' from x^y^y' 

This means bringing the representative matrix and component matrix into the form 

/ S[j 0 0 0 \ 

0 Sj a a 

0 a s 

0 a s J 


To decouple x' from x, y, y', calculate hrst S' with x decoupled from x', y, y' (see above). 
Then S" = has x' decoupled (Ci is symplectic). 

The decoupling is preserved under the transformation group R = exp + 7 ^) 0 

+ 1 (/^3 + /^ 2 ) Xi + 01 x2 + /^sXs] (free parameters </>,Xi,X 2 ,X 3 )- 


( Sn 

0 

CO 

Sl4 \ 

0 

S22 

0 

0 

S31 

0 

S33 

S34 

V S41 

0 

S43 

S44 j 


(R,G) = -1.43, (R,B) = -4.41, {G,B) = 2.18 



(R,G) = -1.43, (R,B) = -4.41, (G,B) = 2.18 



Figure 25: 

Beam with x' decoupled. 
Coefficient matrix = 


/ 4.4 

0 

0 

0 \ 

0 

2.6 

0.2 

0.8 

0 

0.2 

2.3 

0.4 

V 0 

0.8 

0.4 

2.8 / 


Decoupling y from x,x',y' 


This means bringing the representative matrix and component matrix into the form 

0 0 \ 


S = 


/ Sii 

S12 

0 

Si 4 

\ 

S21 

S22 

0 

S24 


0 

0 

S33 

0 


V S41 

S42 

0 

S44 

/ 


(S^) = 


/ 

0 
0 
0 


0 

si 

a 

s 


y2 

^2 


S 

a 

y3 

^3 


To decouple y from x, x', y', calculate first S' with x decoupled from x', y, y' (see above). 
Then S" = Cl has y decoupled (C 2 symplectic). 

The decoupling is preserved under the transformation group R = exp[pCi-7')'^ 
+ ^ (/^3 ~ (^2) Xl + P2X2 + f^lxs] (free parameters 0 , Xi, X2, Xs)- 
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(R,G) = -1.43, (R,B) = -4.41, {G,B) = 2.18 



(R,G) = -1.43, (R,B) = -4.41, (G,B) = 2.18 



Figure 26: 

Beam with y decoupled 
Coefficient matrix = 


/ 4.4 

0 

0 

0 \ 

0 

2.6 

0.2 

0.8 

0 

0.2 

2.3 

0.4 

V 0 

0.8 

0.4 

2.8 / 


Decoupling y' from x,x',y 


This means bringing the representative matrix and component matrix into the form 

0 0 \ 


E = 


( Sn 

S12 

Sl3 

0 

\ 

S21 

S22 

S23 

0 


S31 

S32 

S33 

0 


V 0 

0 

0 

S44 

/ 


(S^) = 


/ 

0 
0 
0 


0 

si 

s 

a 


y2 


a 

a 

y3 


To decouple y' from x, x', y, calculate first S' with x decoupled from x', y, y' (see above). 
Then S" = (^gS'^g has y' decoupled {Cs is symplectic). 

The decoupling is preserved under the transformation group R = exp[pCi-7')'^ 
+ ^ (/^3 ~ Xi + 01x2 + 0lx^ (free parameters 0, Xi, X 2 , X 3 )- 


(R,G) = -1.43, (R,B) = -4.41, (G,B) = 2.18 



(R,G) = -1.43, (R,B) = -4.41, (G,B) = 2.18 



Figure 27: 

Beam with y' decoupled. 
Coefficient matrix = 


/ 4.4 

0 

0 

0 \ 

0 

2.6 

0.2 

0.8 

0 

0.2 

2.3 

0.4 

V 0 

0.8 

0.4 

2.8 / 
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A Calculating the normal decomposition with Octave 


yX Octave test program (should run on Matlab, too) 

yXo Create random beam matrix and find normal decomposition 

gam=[0 1 0 0 ;-l 0 0 0;0 0 0 1;0 0-10]; 
bet=abs(gam);Q=(eye(4)+i*bet)/sqrt(2); 
sig=randn(4);sig=sig+sig'; 

y. Step 1 
bo=(0==0); 
while bo 
sig=sig+eye(4); 

[E,L]=eig(sig*gam); 
bo=(norm(real(L))>lE-12); 
end 
sig 

dispC’step 1: eigenvalue decomposition’); disp(’ ’); 
diagL=diag(L). ’ 

E=E 


y Step 2 

dispC’step 2: eigenvectors ordered’);disp(’ ’); 
L=bet*L*bet; 
diagL=diag(L).’ 

E=E*bet; 

y Test: should be 0 
testnil=norm(E*L/E-sig*gam) 
disp(’ ’); 

y Step 3 

dispC’step 3: eigenvectors symplectic’);disp(’ ’); 
h=E.’*gam*E; 

E=E/diag(sqrt([-h(2,1),h(l,2),-h(4,3),h(3,4)])) 
yiest: sollte gam sein 
testgam=round(100*E.’*gam*E+lE-12)/100 

y Step 4 

dispC’step 4: calculate N’); disp(’ ’); 

N=E*Q 

y Test: should be gam 
testgam=round(100*N*gam*N.’+1E-12)/100 
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B Supplements to the Clifford algebra C/3^i(M) 


B.l How to represent the units of 


The 16 ’’real Dirac matrices” 1, 7 ^, /3^ representing the nnits of Cls^i{M.) can be 

constrncted by forming tensor prodncts of the 4 ’’real Panli matrices” l,/ 3 ]^,/ 32,7 
representing the nnits of C/ 2 ,o(®)- The tensor product of two 2x2 matrices is 


ail 

ai2 

«21 

0-22 


bn bi2 

b21 b22 


an 

611 

^21 

bi2 

b22 

012 

'bn 

bi2 

b22 



011611 

011612 

012611 

012612 



^21 



011621 

011622 

012621 

012622 


'bn 

b2i 

bi2 

b22 


'bn 

&21 

bi2 

b22 



021611 

021612 

022611 

022612 

021 

O22 



.021621 

021622 

022621 

022622. 


The matrix multiplication law is (a 0 6 )(c < 8 ) d) = (ac 0 bd). 

There are 24 ways (not counting sign change) to represent the generators 7 ^, /3^, f3l, 
I3\. The representation of the remaining 12 real Dirac matrices then follows from the 
multiplication table in section 5.1. A possible (and natural-looking) representation is 


1 


7^ 


V 



.c 


p\ 






2 x 2 

10 1 

10 7 

7 0/3i 

7 0/3 

0 

-/3i 0 7 

/3i 0 1 

^ 2 ® (^2 

-(32® (3 

2 x 2 

-701 

-707 

10 /3i 

10/3 


-/32 0 7 

/32 0 1 

-/3i 0 (32 

/3i 0/3 


In ref.s n. 0. the units are named 7 „, n = 0... 15. The correspondence is as follows: 


CO.H 

7o 7i 

72 

73 

74 

75 

76 

77 

78 

79 

7io 

In 

712 

713 

714 715 

here 

7^ /3^ 

01 

-f^l 

1 

CO CO 

-/3? 

/3^ 

-c" 

-c^ 



- 0 , 

-01 

02 

1 

A 

CO 


B.2 Calculating the effect of an elementary transformation 

The calculation of the transformed beam matrix S' = RSRt under an elementary 
transformation requires the multiplication table of the algebra units given earlier, but 
also the following table of ’’two-sided products” of the form ABA. In this table, 
e = ( 61 , 62 , 63 )''' and (/^,/^,/^) are unit-vectors. 



■ 1 

1 


Ck 

/3r 


1...1 

1 


Ck 

Pk 

PY- 

■ ■Yr 

-1 

_ 2 

— Ck 

f3k - 2 r(r/3D 

6 - C 

. . .C • 6 

-1 


Ck - 26 fce ■ C 

(3^ - 2e,{e,(3T) 

e^r(3■■■■ 

(3UY 

1 

7 2/'"(/''7”) 

-Cfc + 26 fce-C 

PC + Are,{rei(3() 
-2fYrPk) - 2ek{eiPC) 
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By way of example, these are the calculation steps for the 7 ^-rotation acting on S; 


S' = R^SRT 


= exp (7V) (lS° + S"*) exp (-7^ 

= ISq + (cos 0 + 7 ^ sin 0 ) ■ S"^ j (cos 0 — 7 ^ sin (f) 

= ISg + cos^ S™ + cos 0 sin 0 7\/9™' 


•S' 


sm 


7I/3-7I) ■ S 


= ISg + COS^ < 


sm 


. S^ + ^ 2 . ^2 + ^ 3 . gsj + 2 cos 0 sin 0 S^ + ^ 2 . gs 


= lS^ + /3i-Si+/3 


cos {2(f)) S^ + sin ( 20 ) S^ 


+ (3^ 


— sin ( 20 ) S^ + cos ( 20 ) S^ 


The calculation of the other three elementary transformations runs similar. For the two 
boosts, sin and cos are replaced by sinh and cosh, and the commutator bracket [, ] is 
replaced by the anti-commutator bracket { , }. 


B.3 Comparison of four diagonalization recipes 

In ref.s m. 0. Baumgarten proposes two diagonalization recipes , which are different 
from each other, and from the two recipes given in section 5.6-7. The four recipes may 
be summarized as follows: 


step 0 ^ 

ref. P, p. 17 

ref. | 2 ], p. 6 

section 5.6-7 

section 5.7 

1 

C2,C3^S3|| 1 

7^ ^ E^T E^ 

^2 ^ s^T E^ 

^2 ^ s2t E^ 

2 

^ E3 = 0 

C25 C 3 1 1 

^3 ^ E^T E^ 

^ E^T E^ 

3 

C^E2 = 0 ,Ei|| 1 

(3l S^T 

C2,C3^Sl|| 1 

7I ^ S2T E3 

4 

(3l ^ E'^ II A; 

7^ ^ E^T E^ 

7I ^ E^T E3 

C ^ E^ll A: 

5 


Cl ^ S^ll A; 

Cl ^ S^ll A: 



Here, = 2 A — EqE^ j, and E^ || k is shorthand for E^ || 1-axis, || 2-axis, 
E^ II 3-axis. 


The reader is challenged to try to understand the geometric ideas behind the 4 recipes, 
and make his own choice - or make up his own recipe. 
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